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For a relativistic charged particle moving in a constant electromagnetic field, its velocity 4- vector 
has been well studied. However, despite the fact that both the electromagnetic field and the equa- 
tions of motion are purely real, the resulting 4- velocity is seemingly due to a complex electromagnetic 
field. This work shows that this is not due to some complex formalism used (such as Clifford algebra) 
but is intrinsically due to the fact that the o(3, 1) Lie algebra of the Lorentz group is equivalent to 
two commuting complex su{2) algebras. Expressing the complex su{2) generators in terms of the 
boost and rotation operators then naturally introduces a complex electromagnetic field. This work 
solves the equation of motion not as a matrix equation, but as an operator evolution equation in 
QQ ' terms of the generators of the Lorentz group. The factorization of the real evolution operator into 

' two commuting complex evolution operators then directly gives the time evolution of the velocity 

, 4-vector without any reference to an intermediate field. 
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1 The equations of motion of a negatively charged particle q = —e moving in a constant electromagnetic field F^^ is 
[ given by 

■ du,, e 

9^: ^ = (^-^^ 
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^-H • where = (xo,r), xq = ct, = (iio,u) and r is the proper time. Once (jl.ip is known, the integration of (fTT 
K*" \ is straightforward. This is a well studied problem, with many published solutionsi'^i^i'^i^i^iii^. Since (|l.ip is just 
■ a matrix equation, it can be directly integrated, as done by Taub^, Hellwig^ and Hyman^. However, the algebraic 
manipulations used were purely formal and gave no insight on why the solution is the way it is. Bacry, Combe and 
Richard* have offered a group theoretic analysis of the problem, but the explicit form of the 4-velocity solution was 
not given. More recently, Hestenes^, Salingaros^, Baylis and Jones^ and Zeni and Rodrigues^ have applied Clifford 
algebra techniques to solve this problem. Salingaros has given a detailed comparison with earlier results of Taub and 
Hestenes, correcting some discrepancies. However, his explicit solution is only for a particle starting from rest. Baylis 
, and Jones' solution for an arbitrary initial 4-velocity is very abbreviated. Zeni and Rodrigues' solution is similar in 
form to the solution presented here, but the derivations are completely different. 

Salingaros' work highlighted the fact that, despite the equation of motion (|l.ip and the electromagnetic field are 
both real, the resulting 4-velocity is seemingly due to a complex electromagnetic field. Since complex fields are a 
natural part of Clifford algebra^iiiSii^, it is unclear whether the solution merely reflect the formalism used or that a 
5^ : complex electromagnetic field is an intrinsic part of the solution. 

In this work, we solve p.ip by an entirely elementary method without invoking any advance formalism such as 
Clifford algebra. Instead of solving (jl.ip as a matrix equation, we solve it via an evolution operator, as it is done in 
the Poisson bracket formulation of symplectic integrators^ "''^^!^^'^'^ '^^. Just like the quantum evolution operator, which 
acts on any quantum state and evolve it forward in time, the classical evolution operator acts on any dynamincal 
variable and evolve it forward in time. The method of classical evolution operator is described in Section II. For 
a relativistic particle in a constant electromagnetic field, the corresponding evolution operator is the exponential of 
generators of the Lorentz group. Since the o(3, 1) Lie algebra of the Lorentz group generators is equivalent to two 
commuting complex su{2) algebras, the evolution operator can be factored into two complex evolution operators. In 
each complex evolution operator, reexpressing the complex su[2) generators as boost and rotation operators then 
naturally introduces a complex electromagnetic field. Thus the complex electromagnetic field is an intrinsic part of 
the solution due to the structure of the Lorentz group, and is independent of any Clifford algebra formalism. This is 
shown in Section III. The explicit velocity solution, for an arbitrary initial 4-velocity vector, is given in Section IV. 
Some conclusions and applications are indicated in Section V. 
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II. CLASSICAL EVOLUTION OPERATOR 

Ea. (|l.ip can be written out explicitly as 



du e ,^ ^ . ,„ ^ , 

— = — (Bxu-Euo) (2.1) 

dr mc 

dun e , ^ , , 

-7^ = — (-E-u). (2.2) 
clT mc 



The evolution of any dynamical variable W{x^,Ufj_) (including a;^ and themselves), is given by 

dW _ dW dXf, _^ dW du^ 
dr dx^ dr du^ dr 



(2.3) 



d e „ 9 



u,F,^—]W^{T + V)W. (2.4) 



oXf^ mc ou^ 

If dXf^/dr and du^/dr are derivable from a Hamiltonian, then the RHS of (|2.3p is just the Poisson bracket. However, 

as realized in Ref.^®, as long as one has the equations of motion, one can define an evolution operator regardless 
whether one has a canonical Hamiltonian. Eq. (|2.4|) has the operator solution, 

W^(x^(T),u^(r)) =e^(^+^)M^(a;^,u^). (2.5) 

In particular, because Ff^^ is constant, (T + V^)"Mp — V^u^, the 4- velocity is evolved simply by 

«^(T)=e"(^+^)«^ = e^V. (2.6) 

From the explicit forms (|2.1|) and (j2.2p . one sees that 



tV = — Bxu . — -E-u- uoE- — 

mc V ou ouq ou 

= e(B-J + E-K) (2.7) 

where ^ = er/mc is the scaled proper time, and 

d d d 

J = ux— K = -u- uo^ (2.8) 

ou ouq ou 

are the three rotation and three boost generators of the Lorcntz group. If E = 0, then 

uo(0 = exp (CB ■J)uQ = UQ (2.9) 



exp ( ^B{B X u) ~ ) u (2.10) 



u(0 = exp (CB • J) u 

11^ . 

du 

where B = H/B is the unit vector of the magnetic field. Using the above form (|2.10p . as shown in RefJ^, the 
expansion of the exponential then directly gives the finite rotation 

u(^) = u + sin(^S)(B X u) + (1 - cos(^S))B x (B x u). (2.11) 

Similarly, if B = 0, then the exponential operator gives 

uo(0 = exp(^E • K)mo 

= cosh(^£:)wo - sinh(^£;)E • u (2.12) 
u(0 = cxp(CE-K)u 

= u + E(^(cosh(^£;)-l)E-u-sinh(^£;)uo). (2.13) 



3 



Eqs. (|2.12p - (|2.13p correspond to a general Lorentz boost in the direction oi F, — E/E preserving the Minkowski norm 
Mg(^) — |u(^)p — Uq — I up. The reason for the choice oi q — —e is that the rotation is then right-handed and the 
Lorentz boost is the standard transformation rather than its inverse. 

If both B and E arc non- vanishing, then the general solution for the velocity 4-vector is given by 

"M(e)=e«(^-'+=^-^)^.,. (2.14) 

From the study of symplectic integratoraiSiiiiiSiiHiii^ such an exponential of two operators can be approximated to 
any order in ^ via the product decomposition 

gaB J+E.K) ^ -Q g6.5B.Jge,5E.K (2.15) 
i 

with suitable coefficients {bi,ei}. Since the effect of each exponential is known (|2.9p - (l2.13[) . the the general product 
can be computed in sequence. The two elementary second order approximations are: 

g{(B.J+E-K) ^ g(C/2)B.Jg5E.Kg(C/2)B.J ^ q(^^3-^ 

^ e(€/2)E.Kg5B.Jg(C/2)E.K ^ 0(^3)^ (2.16) 

This approach is purely real, no complex quantity enters anywhere. However, since we are interested only in the exact 
evaluation of p.l4p . higher order approximations of the form (|2.15p will not be considered here. 

III. COMPLEX DECOMPOSITION 

The Lorentz group generators (|2.8p obey the well known commutator relations: 

[Ji,Jj] — -SijkJk [JiiKj] = —SijkKk [Ki,Kj] = SijkJk- (3.1) 

This implies that 

[B • J,E -K] = (-B X E) • K. (3.2) 
Thus if B X E = 0, then B • J and E • K commute and the solution is given by 

"m(0 =e«^-'e«=^''"M (3-3) 

This is the common starting point of many solutions^'^. The remaining task is then to transform to such a Lorentz 
frame in which B x E = 0, apply the solution p.3p and transform back^>^. However, there is a much more direct 
way of decomposing B • J + E • K into two commuting operators based on the fact that the Lie algebra of o(3, 1) is 
equivalent to that of su{2) x su{2). The complex operators 

M=i(J + iK) and M* = i(J-iK) (3.4) 

have commutators, 

[M„ M,] = -EyfeMfe [m;,m;] = -£yfeii/4* [m„m;] = o. (3.5) 

Thus we have an exact, but complex decomposition 

gC(B-J+E-K) ^ ^^(B+iE)-M*gJ(B-iE)-M_ (g g^ 

In p.6p . one notices that 

gC(B-iE)-M ^ gC(B-jE)-(J+»K)/2 (3 j-j 

corresponding to having a complex electric field F/2 = (E + iB)/2 and a complex magnetic field —iF/2 = (B — iE)/2. 
Since —iF x F = 0, one can further decompose, 

gC(B-jE)-(J+jK)/2 ^ g{(F/2)-Kg4(-jF/2).J^ (3 3^ 
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Similarly, 

g5(B+iE).(J-jK)/2 ^ ^C(-jF/2)*-Jg^(F/2)*-K^ (^3 9^ 

Thus the complex decomposition of the purely real evolution operator exp[^(B • J + E • K)] naturally introduces the 
complex electromagnetic field F = E + iB. Furthermore, the evolution operator consists of four pieces, best evaluated 
in pair as in (|3.8p and (|3.9p . Note that 

g€(F/2)*.KgC(F/2).K_^gCE.K (g^^p) 

and the evaluation of ()3.8|) and (|3.9p cannot be further simplified. The complex decomposition has been used by 
Fredstedi^ to obtain a spinor representation of the Lorentz group, but he did not solve the charged particle problem 
as it is done here. 

IV. EXPLICIT SOLUTION 

In order to evaluate (|3.8p and (|3.9p using our previous results p.9p - (l2.13p . it is only necessary to define the norm 
of the complex vector F = E + iB as 



F = VF • F = v/|E|2 - |B|2 + i2(E • B) (4.1) 
and the corresponding complex unit vector as F = F/F. For z = x + iy = \z\e^^ , ^/z = u + iv, where 

u = \z\^^^cos{e/2) = ^|z|(l + cos6l)/2 = y/{\z\+x)/2 

V = |z|i/2sin(6i/2) = -cos6i)/2 = V(|z| -a;)/2. (4.2) 
For the ease of later comparison, we will follow Salingaros' notation and define 



Ki |E|" - IB]-' K2 = 2(E-B) K^Jni + ni (4.3) 



Thus 

F = E' + iB' and F*F ^ {E' f + {B' f = k (4.5) 
Since two exponential operators in (|3.8p commutes, their effects simply add: 

^ = g«(F/2) Kg4(-»F/2)-J ^ Wo ^ 

resulting in 

Uo = cosh(^F/2)iio - sinh(^i^/2)F • u 

u' = u + F(^(cosh(CF/2)- l)F-u-sinh(^i^/2)Mo) 

+ sin(-i$i^/2)(F X u) + (1 - cos(-iCi^/2))F x (F x u) (4.7) 

The latter now greatly simplifies to 

u' = cosh(^F/2)u - sinh(^F/2)(MoF + iF x u). (4.8) 

One can check that the cross-product term above is essential for preserving the Minkowski norm (uq)^ — (u')^ = u^ — v?. 
To further simplify the notation, let's denote 

c = cosh(^i^/2) and s = sinh(CF/2) (4.9) 
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then, finally, 



where 



MO 

u(0 



g5(-»F/2)*-Jg4(F/2)*.K I Wq 



uoiO = c*u'f^ - s*F* ■ u' 
u(0 



(c*c + s*sF* • F)uo ~ (2Re{c*sF) - is*sF* x f) • u 
c*u' ~ s*{u'„F* -iF* X u') 
= (c*c - s*sF* • F)u - (^2Re(c*sF) + is*sF* x f) uq 
+2Im(c*sF) X u + s*s(F*F + FF*) • u. 
All terms can now be easily computed: 



c c : 



cosh(^£;') cos(^B') 



s s 



cosh(^£;') cos(^B') 



(4.10) 



(4.11) 



(4.12) 



sinh(g^O ^ . sinj^B') 



F = e + ib 



e = - (£;'E + B'B) b = - (£;'B - B'F) 

K K 



F* F 



F* X F = 2i- 



E X B 



F*F + FF* 



EE + BB 



E B 

2Re(c*sF) = -[E' sinhiCE') + B' sini^B')] + —[B' smh{^E') - E' sm{^B')] 

K K 

2Im(c*sF) = -[£;'sm(^B') - B' smh{^E')] + — [B'sm(^B') + sinh(^£;')]. 

K K 

Putting everything together then yields 

'cosh(^£;') cos(CS') E^ + B"^ rcosh{(E') coa{^B') 



MO - (- 



Z Z J / 

-[£;'sinh(C£;') +B'sin(^B')] + —[B'smh{^E') - E' sin(^B')] 

K K 



E X B 



[cosh(C£;') -cos(CB')]) -u 



u(0 



/cosh(C£;') cos(fB') ,^^2 rcosh(f£:') cos{^B') 

\ 2 ' 2 K 2 2 



u 



-[£;' sinh(C£;') +B'sin(^B')] + — [S' sinh(^£;') - E' sin(^B')] 



cosh{^E') ~ cos{^B')])uo 



E X B 



+ (-[£;' sin(CB') -^'sinh(CS')] + —[B' sm{^B') + E' smh{^E')]J x u 
+ [cosh(ei^O-cos(OT^(^-")+^(^-"). 



(4.13) 



(4.14) 



The correctness of (|4.13p and (|4.14p can be checked in four ways. First and second, for E = and B = 0, (|4.13p - (|4.14p 
reproduces ((2:9)) - ((2Tl|) and ((2T^ - ((2l3l) respectively. Thir d, for Ex B = 0, fiT^ and ((4T4)) are identical to ([SlSjl . 
Fourth, if initially u = and uq ~ 1 {c = 1), then ()4.13p and ()4.14p arc in complete agreement with Salingaros's 
result^, his Eq.(28). (Recalls that we take q = — e, ^ = eT/m here is the negative of Salingaros' ^ = qT/m.) Thus 
every term in (I4.13P and (|4.14p has been verified in at least one special case. 
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V. 



CONCLUSIONS 



In this work, we have solved the problem of relativisitc motion in a constant electromagnetic field by the method 
of classical evolution operator. The distinct advantage of this approach is that the resulting evolution operator is just 
the finite group generator of the Lorentz group and the fundamental group structure of the Lorentz group can be 
fully exploited for its evaluation. The resulting complex decomposition explains why the solution, despite being real, 
is fundamentally due to a complex electromagnetic field. 

In contrast to other methods, which can only solve the problem when the electromagnetical field is constant, our 
operator approach can be extended to solve the case when the electromagnetic field is no longer spatially uniform. In 
this general case, the full evolution operator (|2.5|) can still be approximated to any order via 



for suitable sets of coefficients {ti,Vi}. Since the effect of both e'^^ and e^^ are known, any dynamical variable W 
can be evolve forward in t by sequentially updating and via 



(5.1) 



u^[viT) = e 



ViTV 



(5.2) 



and 



(5.3) 



The detail of solving relativistic motion in a non-uniform field will be given in a future study. 
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